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Abstract. - We discuss the response of a quantum system to a time-dependent perturbation 
with spectrum <&(u). This is characterised by a rate constant D describing the diffusion of 
occupation probability between levels. We calculate the transition rates by first order per- 
turbation theory, so that multiplying <f>(w) by a constant A changes the diffusion constant to 
XD. However, we discuss circumstances where this linearity does not extend to the function 
space of intensities, so that if intensities <J?i(w) yield diffusion constants Di, then the intensity 
does not result in a diffusion constant J^-Dj. This 'semilinear' response can occur 
in the absorption of radiation by small metal particles. 



1. Introduction. We describe a previously unremarked phenomenon concerning the re- 
sponse of a quantum system to a time-dependent perturbation. The effect is significant if the 
characteristic frequency scale of the perturbation loq obeys ghuo < 1, where g is the density of 
states. Our analysis is an extension of linear-response theory, in that it too relies on first-order 
perturbation theory, in our case to derive rate constants for transitions between levels. The 
rate constants are used in a master equation, which is then analysed non-perturbatively. The 
response obtained is always linear in the intensity of the driving perturbation, but in some 
circumstances the response is not a linear functional of its spectrum. This point will stated 
more precisely below. Conventional linear response always describes the initial response of a 
system, whereas our theory also considers how the response may differ after an initial tran- 
sient. The data plotted in figure ^ (which is explained later) show that the predictions can 
differ by orders of magnitude. 

We characterise the response of the system by the rate E = dE/dt at which its energy 
is increased by the action of the perturbation. This is directly related to experimentally 
observable quantities, such as the absorption of radiation by small metallic particles in an 
electromagnetic field. The expectation from conventional linear-response theory (see for ex- 
ample [1]) is that the rate of absorption is a linear functional of the spectral intensity 
of the radiation: 

/>oc 

E = duj a(u)$(u) (1) 
Jo 

where a{uj) is a frequency-dependent absorption coefficient. This expression satisfies two 
requirements for linearity: if an intensity function $i(w) results in a response E^, then Q 
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implies that for some constant A 



$(w) i-> A$(w) 



E^ XE 



(2) 
(3) 



In this letter we introduce a new form of linear-response theory, which satisfies (|5J) but not 
@, and which we therefore term semilinear response theory. In the limiting case gftioo <C 1 
we show that the absorption rate is approximated by 



where will be specified later. Equation clearly satisfies condition J3J), but not |J3}. This 
is a consequence of the fact that Q is a weighted harmonic average of the spectral intensity 
function, $(w). Reference [2] discusses a model for a quantum dot coupled to a conducting 
ring, where the DC conductance is also obtained by a harmonic averaging procedure. In 
that case, however, harmonic averaging is relevant because of the specific structure of the 
Hamiltonian of that system. By contrast, the results described here relating the AC response 
to <E>(w) are applicable to generic systems. 

Our approach is based upon an observation about the response of quantum systems to low- 
frequency perturbations. We discuss a system which absorbs a finite amount of energy, E , 
and consider taking the limit as the characteristic frequency ujq of the perturbation approaches 
zero. The number of quanta that the system absorbs, E /huj , diverges as oj — > 0. In order 
to understand the response of a quantum system to low-frequency perturbations, we must 
therefore consider multiple excitations. We describe the excitation of the system by a master 
equation, describing the probability p n (t) that the system is in state with level number n 
at time t. If this master equation is treated in perturbation theory, we recover conventional 
linear-response theory. However, our non-perturbative treatment yields distinctive differences 
from the usual linear-response results. 

Our results are quite generally applicable, but it may be helpful to bear in mind a specific 
example, namely a single electron trapped inside an irregularly shaped enclosure, subjected 
to fluctuating electric fields. This is a simplified model for the absorption of electromagnetic 
radiation by small metallic particles. In a classic paper, Gorkov and Eliashberg [3] predicted 
a quantum size effect, where the absorption of radiation would show distinctive structures 
for frequencies close to the frequency (gfr) -1 (here g is the density of states of single electron 
excitations at the Fermi energy). The energy levels were assumed to have the same statistical 
properties as random matrices, and the absorption was calculated using random matrix models 
introduced by Dyson [4], discussed in [5]. We start from the same model and arrive at 
very different conclusions. We remark that despite intensive investigation, there is no clear 
experimental evidence for the validity of the theory proposed in [3] . 

2. The Hamiltonian. We denote the Hamiltonian in the absence of external fields by Ho- 
The perturbation is described by a set of N p operators Hj multiplied by time-dependent fields 



In the case where the theory is applied to very small metal particles in an electromagnetic 
field, TCq is the Hamiltonian for quasiparticle excitations, and the TLj are operators repre- 
senting coupling of quasiparticles to the -/V p = 3 components of the electric field, Xj(t). The 
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Fig. 1 - (a) The network of transitions between states is analogous to a random resistor network, 
(b) Rate of energy absorption divided by difference in the occupation probabilities of the lowest 
and the highest states, showing crossover from an initial transient to a steady state. The numerical 
simulations of li lt employ random-matrix eigenvalues with /3 = 1, cr — 1, g = 1, iV p = 3, and 
$H = e 2 exp(-HMOM), for e 2 = 1 and (Jo = 2 (o), 1 (>), 0.5 (0), 0.3 (□), and 0.25 (<). The 
results were averaged over 200 realisations of TLj. The asymptotes are: left, linear-response for the 
initial transient feq. J21l l. right, steady state determined by network model feq. H4l l. (c) Energy- 
diffusion constant De vs. ojq. Symbols are long-time asymptotes from figure 03. The solid line 
is the semilinear response obtained by solving 1141 . averaging over 20 realisations. Also shown are 
linear-response approximation fea. l21H l. asymptote to semilinear-response theory [eg . 1221 1. 



operators Ti.j are not simple dipole operators, because they must take account of screening of 
the externally applied perturbation by polarisation charges [6,7]. 

The fields are not monochromatic, and their components Xj(t) are random functions of 
time, satisfying 

(X J (i))=0 and (X i (t)X J (t / ))=%^(t-i / ) (6) 
(angular brackets denote averages). The fields have a spectral intensity &(u>), defined by 

/oo 
dr 0(r) exp(iwr) . (7) 
-oo 

In the numerical examples below, we take = e 2 exp(— \u>\ / coq) / ojq where Uq, e are con- 

stants. An exponential dependence with loq = k^T /h is a natural choice if the system is 
excited thermally. 

3. Master equation. First consider the rate for transitions between eigenstates of Ho- 
Expand the solution \ift(t)) of the Schrodinger equation ihd\ip)/dt = 7i(t)\ip) 

l^(*)> = S a " ( ^ exp(-iE n t/K) \<p n ) (8) 
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where Ti.o\(p n ) = E n \ip n ) and the energies are ordered according to the index n. The amplitudes 
a n (t) satisfy: 

-i Np 
a n = -z- ^ exp(iw nm t) H^-Xj W a m (9) 
m j=l 

where W^m = (^nl'Wjl^m) and w nOT = — E rn )/h. Solving perturbatively for the initial 
condition a n (Q) = 5 nm one obtains an expression for a n (t) and hence for the probability 
Pn(t) = {\a n (t)\ 2 ) to be in the n th state. For n ^ m we have p n (i) = r nm t + 0(t 2 ), where the 
rate constants are given by a version of Fermi's golden rule 

1 Np 

r nm = ¥ $( Wnm ) £|w«| 2 . (io) 

The expression for p„ is valid for times sufficiently short that 1 — p m <C 1, but large enough 
that u>ot ^> 1. These conditions are compatible for a sufficiently small e. 

The rate constants can be used to write the master equation for the occupation probabilities 
Pn(t): 

-TjT = ^2^nm{Pm ~ Pn) ■ (H) 

m 

This master-equation ignores interference effects, which average away on long timescales. 

Our model © describes excitation of a quantum system due to a perturbation. In 
practical applications, the system may also be subject to relaxation effects. The master- 
equation model can then be augmented with terms representing relaxation processes. These 
additional terms could represent the transfer of energy from electronic excitations into phonons 
or photons. Electron-electron interactions could also be included, although these represent 
re-arrangement of energy within the electronic system rather than relaxation. The lifetime 
for an electron to emit photons or phonons diverges as the energy of excitation of the electron 
approaches zero. For systems excited by low-frequency fields, electrons may therefore be 
excited by many quanta before their relaxation rate is significant. Thus our approach is in 
contrast with conventional linear response theory [10] which implicitly assumes that multiple 
excitations are not relevant. 

In the case of absorption of electromagnetic radiation by small conducting particles, ref- 
erence [3] shows that phonons do not cause relaxation at low frequencies, so that emission 
of photons is the dominant relaxation mechanism. In this case it is easy to see that there is 
multiple excitation when the intensity of radiation at frequency ujo is large compared to the 
intensity of black-body radiation at temperature T — hujo/kB- This condition is easily sat- 
isfied at the microwave or far-infrared frequencies which are relevant to experimental studies 
on the effect discussed in [3]. 

4. Energy diffusion and resistor networks. We are interested in the long-time behaviour 
of the master equation (|llfl . The coarse-grained occupation probabilities obey a continuity 
equation: 

d tP (n,t)+d n J(n,t) = Q (12) 

with probability flux J. We argue below that the coarse-grained occupation probability obeys 
Fick's law, J = —Ddp/dn, so that p(n,t) obeys a diffusion equation 



d t p = d n [Dd n p] 



(13) 
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In order to determine D we make use of the analogy between (jllf) and Kirchoff 's equation for 
a resistor network (illustrated in figure ^^) : nodes n and m are connected by conductances 
G nm . If a current is supplied at node n, the potentials V n satisfy 

4=^G nm (K-y m ). (14) 

m 

The probabilities p n in I jlljl correspond to the potentials V n , the rates T nm to the conductances 
G nm , and in the steady state I n = —dp n /dt = at all the nodes. Coarse graining is achieved 
by considering a truncated network segment of length An > 1 as illustrated in figure QJi, 
with a current J is injected into one end and extracted at the other end. This finite segment 
is described by an equation in the form of l|14l) with /„ = J5 n ,n ~ J3n,n +An> resulting in 
a potential difference AV — V no +& n — V no . We expect that AV/An approaches a limit as 
An increases, implying a 'coarse grained' Fick's law: J « — D(p na+ &„ — p no ) / An, which is 
analogous to Ohm's Law. Thus D j An is the conductance of the segment, and D is obtained 
as D = — JlimAn^oo An/ AV. We remark that Miller and Abrahams [8] introduced random 
resistor models in studies of spatial (as opposed to energy) diffusion in disordered systems. 

In what follows we discuss expressions for D in two limiting cases. When the rate constants 
are negligible for all but nearest-neighbour transitions, the resistance (D/ An)' 1 is the sum 
of resistors in series, leading to 

nj+An s —1 

1 \ " ri-1 1 — An-l 



n— no 



To show that Fick's law holds in the limit considered here, it suffices to show that (|15|) yields a 
finite result for D (an example is given in equation 1221) below) . It is an immediate consequence 
of the nature of the harmonic average that the diffusion constant is significantly reduced by 
the presence of 'bottlenecks', that is links with very low transition rates [9]. Reference [2] 
considers the consequences of this for energy diffusion in a model for the DC response of a 
quantum dot coupled to a conducting ring. 

We turn to the other extreme case, where many transitions have significant weight (not just 
to near neighbours). We may assume that the potential changes linearly along the network 
(p n oc n) and find that 

D = ( 5^(™-") 2 rn,m)- (16) 

m 

where the angular brackets means averaging over n as in equation (|15|) . To derive this, note 
that the contribution of bond of length An = (m — n) to the current is proportional to the 
potential drop and hence to An. Furthermore the number of bonds of length An passing 
through a given section gives a further factor An/2. 

5. Rate of absorption of energy. The expectation value of the energy of the system is 

E{t)=J2Pn(t)E n . (17) 

n 

When many states are excited, the sum may be approximated by an integral. Applying (|13|l . 
the rate at which energy is absorbed is 

E = Jdn E n d„ [Dd n p] =-JdEgD E ^. (18) 
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This equation shows that the rate of energy absorption E is proportional to the energy- 
diffusion coefficient De = D/g 2 , a principle that was introduced in [11]. 

6. An example of semilinear response. Consider the case where ghujQ <C 1. Here the rate 
constants (| 1 0|) decrease very rapidly as the separation in energy increases, and we can neglect 
all of the rate constants T nm other than those describing nearest neighbour coupling. The 
diffusion constant De = D/g 2 is then estimated via equation l)15[l. It is clear that large gaps 
in the spectrum create 'bottlenecks' which slow the diffusion of probability. We define P(S)dS 
as the probability that the normalised spacing between two successive levels (E n +\ — E n )g 
is in the interval [S, S + dS] . We assume that the matrix elements Hum are independent 

Gaussian random variables with variance a 2 and zero mean, independent of the energy levels. 

(i) (1) (i) 

We write Hnm — crXrtm with Gaussian random variables Xnm, each with zero mean and unit 

variance, then substitute (|10fl into l|15fl and find 

dxe~ x2 / 2 



De = 



{ghf 



27r )JV P /2 x 2 



(19) 



Eq. p9|l is an example of semilinear response, in the form of eq.(0}. Assuming that $(w) 
decreases rapidly when u> ^S> ojq, the integral in l|19|l is dominated by the tail of the level 
spacing distribution. Denoting the second term (with the integral over dx) by W, we find 
that W = 1 for N p = 3. For N p = 1, 2 we find W = 0, so that D E = 0, implying that the 
spread of probability is sub-diffusive. 

7. Linear-response theory. Now we contrast l|19|l with conventional linear-response theory. 
We assume that the initial probability p n (0) is a smooth function of the energy of the state. 
We differentiate (|llfl . substitute in JT7j, and expand p m — p n to first order in E n — E m . 
Interchanging the indices n and m and averaging the two expressions for the double sum 
gives: 

E = — \ ^](£ m — E n ) 2 T nm -^ . (20) 

n,m 

We can identify De by comparison with l|18fl . After substitution of i|l(J (I we obtain an expres- 
sion in terms of the two-level correlation function i?2(e) = J2nm(3(E — E n )8(E + e — E m )}/ g 2 
(we use the notation of [5]) 

POO 

D E = N p (j 2 hg duuj 2 R 2 (huj)<S>(uj) . (21) 
Jo 

This a linear functional of the spectral intensity $(w), leading to and expression of the form QJ. 
Equation (|21|l is a version of the 'Kubo formula' of linear-response theory and is equivalent 
to a result obtained by Gorkov and Eliashberg [3]. It is subject to the criticism that it 
only describes the initial response of the system (figure^)): after a short transient, the 
probabilities may cease to be given accurately by a smooth function of the energy, leading to 
a very different rate of absorption, such as that given by l|19|l . Conventional linear-response 
theory [10] implicitly assumes that strong relaxation prevents level-number diffusion from 
exploring the 'bottlenecks'. 

Finally we remark that if the initial probability p n (0) is a smooth of level number n instead 
of E n , a slightly modified form of linear-response theory is obtained [10]. 

8. Random-matrix models and numerical experiments. We now compare the calculation 
of De using (|19|) with conventional linear- response theory 1)2 l|l . assuming random matrix 
models for P(S) and i?2(e)- Although the spectra of complex quantum systems differ, their 
statistical properties are very similar and can be calculated for suitably defined random-matrix 
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ensembles [5], There are three 'universal' ensembles, labelled by an integer index /3 £ {1, 2, 4}. 
For P(S) we use the 'Wigner surmise', P(S) ~ apS 13 exp(— cpS 2 ), with ap and chosen so 
that P(S) is normalised with mean value unity. The diffusion coefficient i|19|) depends on 
the large separations, so that we require accurate information about the values of P(S) for 
large argument: precise information about the large S asymptotics is given in [5], but for the 
moderately large values of S that are probed by our numerical studies, the Wigner surmise 
gives more accurate results. We also require ife(e), to evaluate l|21|). Here it is the behaviour 
for small spacings that is of most interest, where R%{e) — kp(ge)^ + 0(e l3+1 ) with universal 
constants kp. 

We illustrate the theory by comparing the predictions of (|19f) and (|21|1 for a spectral 
intensity $(w) = e 2 exp(— \w\ / Wo) / ojq and (3=1, using the Wigner surmise for P(S). As 
ujq — > 0, equation 119|) predicts that De approaches zero in a non-analytic fashion: 

D E = (e 2 a 2 /2gh) exp{-l/[^(^ ) 2 ]} . (22) 

This is dramatically different from the result of conventional linear response theory l|21|) where, 
for small values of us , we find De ~ CpN p cr 2 e 2 (hg)^ +1 Wq +2 (for some universal constants Cp). 

Figure \^p,c shows numerical results for simulations using random-matrix energy levels, 
with j3 = 1, e = 1, <7 2 = 1, N p = 3, and g = 1. The data in figure are obtained by a 
simulation of with (|17|l . We use GOE [4] random matrices of dimension iV = 4000 and 
an exponential initial distribution p n oc exp(— E n /AE) with gAE = 100. After an initial 
transient both E and (pi(t) — pj\r(t)) decrease, their ratio approaching a limit which (using 
Ijl8|l ) we identify as gDE- These limiting values of De are plotted as symbols in figure^. 
The data for the resistor network (solid line) was obtained by solving Kirchoff's law for a 
network with N — 4000 nodes, using singular-value decomposition. 

9. Summary. We have considered the non-perturbative solution of a master equation 
describing transitions between levels. Its solutions are in general diffusive for large times, with 
a diffusion constant obtained from the conductivity of a random resistor network. When the 
characteristic frequency ujq is small (ghujQ «C 1), only transitions between neighbouring levels 
are significant, analogous to resistors in series. The diffusion constant is then the harmonic 
mean of the rate constants. It is determined by the tail of the level spacing distribution and 
is an example of semilinear response. 
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